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ABSTRACT 

Click  data  from  a  tagged  Mesoplodon  densirostris  was  compared  with  broadband  acoustic  recordings  from 
an  82  hydrophone  wide-baseline  array  located  in  the  Tongue  of  the  Ocean,  Bahamas.  Two  detectors,  an 
Fast  Fourier  Transform  (FFT)  based  detector  and  matched  filter,  were  evaluated  in  white  noise  and  with  the 
acoustic  recordings  from  the  array  for  performance  detecting  M.  densirostris  clicks.  The  matched  filter 
performed  the  best,  allowing  92%  of  the  tagged  animal’s  clicks  to  be  detected  on  at  least  one  hydrophone. 
Time-Difference-of- Arrivals  (TDOAs)  between  the  DTag  and  the  surrounding  hydrophones  were 
computed.  These  TDOAs  were  used  to  compute  a  three-dimensional  hyperbolic  localization  track  of  the 
tagged  animal.  A  maximum  detection  range  of  6500  m  from  the  tagged  animal  to  the  recording 
hydrophone  was  observed.  Offset  aspect  angles  were  determined  from  the  DTag  heading  information  and 
the  bearing  to  the  receiving  hydrophone.  Clicks  within  ±30  degrees  were  detected  at  the  farthest  ranges, 
while  clicks  were  detected  at  all  off-set  angles  at  closer  ranges. 

Cliquez  donnees  d’un  Mesoplodon  densirostris  a  ete  compare  a  large  bande  a  partir  d’un  enregistrement 
sonore  82  hydrophones  large  eventail  de  reference  situe  dans  la  Langue  de  l’ocean,  Bahamas.  Deux 
detecteurs,  une  transformation  rapide  de  Fourier  (FFT)  basee  detecteur  et  filtre  adapte,  ont  ete  evalues  en 
blanc  avec  le  bruit  et  l’acoustique  des  enregistrements  de  l'ensemble  de  la  performance  detecter  M. 
Blainville  clics.  Le  filtre  adapte  est  la  plus  performante,  ce  qui  permet  a  92%  des  animaux  marques  de  clics 
pour  etre  detectes  sur  au  moins  un  des  hydrophones.  Time  -  Difference  de  Arrivees  (TDOAs)  entre  les 
DTag  et  dans  les  environs  des  hydrophones  ont  ete  calcules.  Ces  TDOAs  ont  ete  utilisees  pour  calculer  une 
tridimensionnel  hyperbolique  voie  de  la  localisation  des  animaux  marques.  Une  gamme  de  detection 
maximale  de  6500  m  de  la  marques  a  l'enregistrement  des  animaux  d’hydrophones  a  ete  observee.  Offset 
aspects  les  angles  ont  ete  determinees  a  partir  de  la  rubrique  DTag  information  et  le  roulement  a  la 
reception  d’hydrophones.  Clicks  a  ±  30  degres  ont  ete  detectees  lors  de  la  gamme  la  plus  eloignee,  tandis 
que  les  clics  ont  ete  detectees  dans  tous  les  angles  hors  serie  a  rapprocher  fourchettes. 


1.  INTRODUCTION 

On  October  23,  2006,  a  Woods  Hole  Oceanographic 
Institute  (WHOI)  DTag  was  placed  on  one  individual  in  a 
group  believed  to  consist  of  four  Mesoplodon  densirostris 
on  the  Atlantic  Undersea  Test  and  Evaluation  Center 
(AUTEC)  underwater  tracking  range.  The  tag  remained 
attached  for  approximately  19  hours  over  which  time  seven 
deep  dives  were  recorded.  Whale  vocalizations  were 
simultaneously  monitored  using  82  AUTEC  bottom- 
mounted  hydrophones.  The  hydrophones  are  deployed  at 
depths  of  -2  km  and  are  separated  by  -4  km  baselines.  Two 
detectors  are  evaluated  for  use  in  M.  densirostris  click 
detection:  an  FFT  based  detector  and  matched  filter. 
Detection  performance  is  first  evaluated  in  the  presence  of 


Gaussian  white  noise,  and  then  with  hydrophone  recordings 
corresponding  to  the  tagging  event. 

2.  METHODS 

2.1  DTag  Data  Description 

The  DTag  was  attached  to  a  probable  female  M.  densirostris 
in  a  group  believed  to  consist  of  two  mother-juvenile  pairs 
at  11:37:38  a.m.  (+/-  5  seconds)  local  time  on  October  23, 
2006.  The  tagging  GPS  location  was  24°  30.412’  N,  IT 
35.320’  W.  A  male  M.  densirostris  may  have  been  in  the 
vicinity.  The  DTag  recorded  stereo  audio  at  a  192  kHz 
sampling  rate  with  an  audio  sensitivity  of  -171  dB  re 
1.0/jLiPa.  The  pitch,  roll,  heading,  and  depth  of  the  whale 
were  determined  from  the  accelerometer,  magnetometer, 
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and  pressure  sensors  sampled  at  50  Hz.  The  DTag 
measurements  were  processed  using  the  methods  described 
in  Johnson  and  Tyack  [  1  ]  resulting  in  processed  orientation 
and  depth  data  with  5  Hz  resolution. 

2.2  AUTEC  Hydrophone  Array 

Prior  to  and  for  the  duration  the  tag  was  attached  to  the 
whale,  audio  data  from  the  82  bottom  mounted  hydrophones 
of  the  AUTEC  tracking  range  were  simultaneously  recorded 
digitally  on  multiple  Alesis  HD24  hard  drive  recorders  at  a 
96  kHz  sampling  rate.  Each  recorder  can  accommodate  12- 
channels  of  data  with  the  last  channel  recording  an  IRIG-B 
modulated  time  signal. 

The  hydrophones  are  mounted  4-5  meters  off  the  sea  floor 
with  an  upward,  roughly  hemispherical,  beam  pattern.  There 
are  68  wideband  hydrophones  with  a  usable  bandwidth  from 
50  Hz  to  approximately  45  kHz.  There  are  an  additional  14 
hydrophones  with  a  bandwidth  from  roughly  8  kHz  to  over 
50  kHz  installed  in  two  7  hydrophone  arrays.  Hydrophone 
data  is  digitized  at  a  sampling  rate  of  96  kHz.  This  is  a 
standard  audio  rate  that  allows  for  Nyquist  sampling  of  the 
wideband  hydrophones,  although  a  folding  effect  is 
introduced  for  the  14  higher  frequency  hydrophones.  This 
folding  has  not  been  found  to  have  a  significant  effect  on  the 
determination  of  click  arrival  times. 

2.3  Detection 

M.  densirostris  produces  echolocation  clicks  with  a 
frequency  modulated  upsweep.  The  peak  source  level  of 
similiarly  sized  delphinids  has  been  estimated  at  220  dB  re. 
1  pPa  [2].  Tag  data  from  another  species  of  beaked  whale 
( Ziphius  cavirostris)  indicate  a  pronounced  beam  pattern 
with  a  3  dB  beam  width  of  6°  [3].  Due  to  the  narrow  beam 
width,  determination  of  the  arrival  time  of  a  specific  click 
on  multiple  hydrophones  in  a  widely  spaced  array  such  as  at 
AUTEC  is  a  challenge  both  for  detection  and  data 
association.  The  hyperbolic  localization  technique  used 
requires  a  minimum  of  three  hydrophones  for  a  two 
dimensional  position  to  be  determined.  Improved  detector 
performance  is  critical  in  order  to  maximize  the  probability 
that  a  given  click  will  be  detected  on  enough  hydrophones 
to  produce  a  position.  Accordingly,  two  detection  methods, 
a  FFT  based  detector  and  a  matched  filter,  were  compared 
with  respect  to  detection  performance.  Detection 
performance  was  first  compared  in  Gaussian  white  noise, 
and  then  on  recorded  data  in  the  vicinity  of  the  tagged 
animal. 

DTag  Click  Detector 

Clicks  recorded  on  the  DTag  were  classified  as  belonging  to 
the  tagged  whale  based  on  two  features.  The  attachment  of 
the  tag  to  the  whale  results  in  a  low-frequency  energy 


component  that  is  not  present  in  clicks  from  conspecific 
whales  [3].  Second,  the  angle  of  arrival  for  clicks  from  the 
tagged  whale  is  close  to  zero  between  the  two  hydrophones 
on  the  tag,  while  it  varies  as  the  whale  moves  for  clicks 
from  conspecific s  [4]. 

FFT  based  detector 

A  multi-stage  FFT  based  energy  detector  has  been 
successfully  used  for  detection  of  clicks  from  a  variety  of 
echo-locating  odontocetes,  including  sperm  whales  [5]  and 
beaked  whales  [6].  A  2048  point  FFT  with  a  50%  overlap  is 
used  for  this  analysis.  At  the  96  kHz  sampling  rate  this 
provides  a  frequency  resolution  (per  bin)  of  46.875  Hz  and  a 
time  resolution  (per  FFT)  of  10.67  ms.  Each  bin  of  the  FFT 
is  independently  thresholded  against  an  exponentially 
decaying  time  average  of  the  data  in  that  bin  as  given  in 
Eq.  1: 

n>0:  NVT[n]  =  {\-a)bin[n]  +  aNVT[n-\ ]  Eq.  la 

n=0:  NVT[n\  =  0  Eq. 

lb 

where,  the  parameter  a  has  been  chosen  empirically  to 
provide  a  time  constant  of  0.2  seconds. 

The  binary  output  of  the  thresholding  process  is  combined 
into  a  single  detection  report.  If  any  of  the  1024  bins  have 
passed  threshold,  the  first  stage  declares  a  detection  and 
passes  the  detection  report  on  to  the  next  stage. 

The  output  of  the  first  stage  of  the  detector  is  then  examined 
to  determine  whether  the  event  was  triggered  by  a  beaked 
whale.  Since  clicks  are  broadband  events,  the  detection 
report  may  be  broadly  classified  as  a  click  by  counting  the 
number  of  FFT  bins  which  triggered.  Assuming  a  click 
event  is  declared,  the  frequency  content  of  the  thresholded 
detection  report  is  examined.  A  set  of  five  frequency  bands 
roughly  conforming  to  species  of  interest  (cut  off  by  the 
hydrophone  response)  have  been  selected,  where  beaked 
whales  comprise  band  2  (Table  1). 


Table  1:  Frequency  bands 


Band 

Low  Frequency  (kHz) 

High  Frequency 
(kHz) 

1 

45 

48 

2 

24 

48 

3 

12 

48 

4 

1.5 

18 

5 

0 

1.5 

A  ratio  of  bins  set  to  the  total  number  of  bins  in  each  band  is 
computed.  If  band  two  is  selected,  then  the  detection  is 
tentatively  classified  as  a  beaked  whale.  Due  to  the  fact  that 
many  of  the  bands  overlap,  a  second  check  is  performed  by 


examining  the  number  of  bins  set  out-of-band.  If  this 
exceeds  10%,  then  the  detection  is  reclassified  as  a  dolphin 
as  they  are  more  likely  to  have  significant  spectral  energy 
below  24  kHz. 

Matched  filter  detector 

A  linear  matched  filter  can  be  shown  to  be  the  optimal 
detector  for  known  signals  in  white  gaussian  noise  [7].  A 
high  signal  to  noise  ratio  M.  densirostris  click  extracted 
from  the  data  set  was  used  as  the  match  template.  The  filter 
implements  a  typical  matched  filter  equation: 

N-l 

mf\p\-  Eq.  2 

k= 0 

where,  mf[n]  is  the  output  of  the  matched  filter  at  point  n, 
x[n-k]  is  the  input  sample  starting  with  the  most  recent  at 
k=0  and  working  backward  in  time,  and  m[k]  is  the  match 
template  containing  the  digitized  click.  N  is  the  number  of 
samples  in  the  digitized  click.  The  instantaneous  output  of 
the  filter  is  then  compared  to  an  exponentially  decaying  time 
average  of  the  filter  output  with  a  time  constant  of  0.1 
seconds.  If  the  instantaneous  output  exceeds  the  time 
average  by  a  specified  threshold,  a  detection  is  declared. 

False  alarm  statistics 


Both  the  FFT  detector  and  the  matched  filter  have  been 
implemented  as  constant  false  alarm  rate  (CFAR)  detectors. 
A  direct  comparison  of  detection  performance  between  the 
two  requires  normalizing  the  false  alarm  rates.  False  alarm 
statistics  have  been  computed  in  the  presence  of  white 
Gaussian  noise  using  the  Box-Mueller  pseudorandom  noise 
generation  algorithm  from  the  GNU  Scientific  Library 
(GSL).  This  algorithm  is  fed  from  a  2nd  generation 
RANLUX  algorithm  with  a  period  of  approximately  10171 
cycles.  The  false  alarm  rate  was  then  computed  by  dividing 
the  number  of  false  detections  by  the  total  run  time  for  each 
threshold. 

Two  sets  of  results  were  compiled  for  the  FFT  detector.  The 
first  set  indicates  the  performance  of  the  first  stage  of  the 
detector.  As  can  be  seen,  this  stage  runs  with  a  high  false 
alarm  rate. 


Noise  Variable  Threshold  (Fixed  Threshold  -100) 


Figure  1:  FFT  Detector  first  stage  false  alarm  curve 

The  false  alarm  rate  drops  dramatically  at  the  output  of  the 
second  stage.  The  main  parameter  determining  performance 
is  the  click  threshold  used  to  determine  when  a  sufficient 
number  of  bins  have  detected  to  declare  a  click  event. 


Beaked  Click  False  Alarms 


Figure  2:  FFT  Detector  False  Alarm  Curve  after  filtering  for 
beaked  whales 

The  matched  filter  curve  is  typical  and  indicates  a  false 
alarm  rate  dropping  exponentially  with  an  increase  in  the 
threshold. 


Mesoplodon  template  at  96  kHz 


Threshold 


Figure  3:  Matched  filter  false  alarm  curve 


Matched  filter  detector,  NVT=28.7009 


A  false  alarm  rate  of  lxl 0"3  was  chosen  as  the  test  criterion. 
To  achieve  this  rate,  the  following  thresholds  were  chosen 
based  on  the  false  alarm  curves: 

FFT  Detector:-35.1432 
Matched  Filter:  28.7009 

Probability  of  Detection 

Probability  of  detection  statistics  were  compiled  by  creating 
a  series  of  test  data  sets  consisting  of  a  high  SNR  click 
identified  as  M.  densirostris.  The  click  was  scaled  by  a 
specified  constant  to  achieve  a  desired  signal  level  and 
repeated  827  times  at  a  regular  interval  of  4/second.  White 
Gaussian  Noise  (sigma=l)  was  added  to  the  signal  to  obtain 
the  desired  SNR. 


FFT  detector,  NVT=-35.1432,  tc=0.2 


SNR  (dB) 

Figure  4:  FFT  detector  probability  of  detection 

The  FFT  detector  tops  out  at  approximately  80%  probability 
of  detection.  Visual  examination  of  the  data  indicates  that 
some  clicks  are  not  present  at  the  output  of  the  first  stage  of 
the  detector.  The  performance  deficit  at  high  SNR  is 
therefore  most  likely  linked  to  the  choice  of  time  constant 
for  the  exponential  noise  filter.  A  lower  time  constant  may 
improve  performance  for  regular  clicks,  however  this 
remains  to  be  investigated. 


Figure  5:  Matched  filter  probability  of  detection.  Two  match 
templates  are  plotted.  First,  template  matching  the  click  used  to 
generate  the  dataset.  Second  template  for  high  SNR  click  from  a 
separate  dataset. 

The  matched  filter  provides  the  expected  behavior  when  the 
match  template  exactly  matches  the  signal  present  in  the 
data  set.  This  is  never  the  case  in  practice.  To  estimate  the 
effect  of  using  an  arbitrarily  chosen  high  SNR  template,  a 
click  from  a  completely  separate  data  set  (also  collected  at 
AUTEC)  was  used  as  a  second  match  template.  This  is 
plotted  in  the  rightmost  curve  in  Figure  5.  Using  an 
arbitrarily  chosen  click  degrades  detector  performance  by 
approximately  2-3  dB.  In  either  case  probability  of  detection 
is  at  least  95%  by  0  dB  SNR. 

Conclusion:  ‘Optimal5  Detector 

The  matched  filter  significantly  outperformed  the  FFT 
detector  on  the  test  data  sets.  This  was  true  even  when  the 
click  used  as  the  match  template  was  not  the  same  as  the 
click  used  to  generate  the  data  sets.  The  probability  of 
detection  data  shows  a  performance  gain  of  at  least  25  dB 
for  the  test  cases  studied.  This  suggests  that  the  click 
structure  is  relatively  constant  with  reasonably  low  variance 
between  clicks  from  different  individuals. 

High  SNR  clicks  are  typically  chosen  as  match  templates. 
Due  to  the  narrow  beam  width  emitted  by  the  animal,  it  is 
expected  that  most  high  SNR  clicks  received  will  be 
received  when  the  animal  has  the  ensonified  receiver 
directly  in  the  beam.  The  structure  of  these  clicks  may  not 
be  representative  of  off-axis  clicks.  It  has  been  suggested 
that  the  click  structure  of  beaked  whale  clicks  is  aspect 
angle  dependent.  In  this  case  the  matched  filter  will  be  sub- 
optimal  at  any  aspect  angle  other  than  the  one  at  which  the 
match  template  was  tuned  for.  However,  in  a  widely  spaced 
array  such  as  AUTEC,  it  is  possible  to  enhance  detection 
and  association  significantly  by  improving  detection 
performance  at  hydrophones  which  are  farther  away  from 
the  animal,  but  still  in  the  beam.  In  this  instance  the  matched 
filter  may  be  employed  to  significant  advantage. 


2.5  Data  Association 


Clicks  originating  from  the  tagged  animal  have  been 
identified  on  the  surrounding  bottom  mounted  hydrophones 
by  matching  inter-click  interval  patterns  [5,8].  These 
patterns  have  been  found  to  be  an  effective  means  of 
associating  patterns  of  detections  among  hydrophones  for 
sperm  whales  ( Physeter  macrocephalus).  A  fundamental 
assumption  is  that  each  animal  exhibits  its  own  unique 
pattern  of  clicks.  The  unique  pattern  is  used  as  a  template 
for  a  comb  sieve  that  is  correlated  against  the  beaked  whale 
clicks  detected  on  the  surrounding  hydrophones  [5].  The 
window  with  greatest  number  of  correlations  between  the 
template  and  the  hydrophone  is  assigned  as  the  time- 
difference-of-arrival  (TDOA)  between  the  DTag  and  the 
hydrophone.  After  the  comb  sieve  is  complete,  the 
probability  density  function  of  the  TDOAs  for  each 
hydrophone  is  calculated  in  one  minute  windows.  TDOAs 
that  are  significantly  above  the  noise  level  in  each  window 
are  passed  on  for  use  in  localization  and  considered  valid. 
The  remaining  TDOAs  are  considered  invalid  and  not  used 
further. 

2.6  3D  Hyperbolic  Localization 

Positions  are  computed  from  the  valid  TDOA  sets  using  a 
hyperbolic  multilateration  positioning  algorithm  developed 
by  Vincent  [9].  Two  2500-ft  depth  XBT  profiles  were 
collected  on  23  October  2006.  These  profiles  were 
combined  with  a  standard  deep  water  profile  and  converted 
to  sound  speed  by  AUTEC.  The  sound  speed  profile 
proximate  to  the  tagging  location  was  used  for  calculating 
the  direct  path  effective  sound  velocity  [9].  TDOAs  are 
required  between  at  least  four  hydrophones  and  the  DTag  in 
order  to  compute  a  3-D  position.  Due  to  the  directional 
nature  of  the  clicks,  there  were  very  few  instances  when  an 
individual  click  was  correlated  on  4  hydrophones. 
Therefore,  the  TDOA  for  each  hydrophone  was  interpolated 
using  the  MATLAB  Piecewise  Cubic  Hermite  Interpolating 
Polynomial  function  (The  Mathworks,  2007)  (Figure  6). 
The  four  hydrophones  with  the  greatest  number  of  valid 
time-difference  of  arrivals,  37,  43,  44,  and  50,  were  used  to 
create  a  Time-of- Arrival  matrix.  The  time-of-emission,  x, 
y  and  z  position  of  the  sound  source  were  estimated  using 
the  TO  A  matrix  in  the  hyperbolic  multilateration  algorithm. 


Figure  6:  Interpreted  TDOAs 


DISCUSSION 

All  results  discussed  in  this  paper  are  for  the  first  deep  dive 
recorded  on  the  DTag,  from  approximately  34  to  91  minutes 
after  tagging.  The  whale  began  vocalizing  6  minutes  into 
the  dive  at  567  m  depth,  and  continued  to  vocalize  for 
approximately  36  minutes  between  567  and  1049  m  depth 
(Figure  7). 


Time  since  Tag  On  (seconds) 

Figure  7:  DTag  Dive  1  Depth  Profile 

3.1  Detection  Efficacy 

Detection  efficacy  was  assessed  by  evaluating  how  many  of 
the  clicks  emitted  by  the  tagged  whale  were  successfully 
detected  and  associated  on  the  nearby  hydrophones.  For 
each  method,  hydrophones  36,  37,  38,  42,  43,  44,  49,  and  50 
were  processed  through  each  detector  and  the  beaked  whale 
filter.  These  hydrophones  were  chosen  by  visually 
evaluating  thresholded  spectrogram  data  over  the  entire 
range  for  the  presence  of  beaked  whale  clicks  for  the 
duration  of  the  tagging  event.  Whale  click  times  from  the 
DTag  were  used  as  the  template  to  search  for  correlation 
with  the  resulting  Time-of- Arrivals  produced  using  the  FFT 
and  matched  filter  detectors.  Using  the  association 


algorithm  detailed  in  Section  2.5,  TDOAs  were  calculated 
between  each  hydrophone  and  the  DTag. 

During  the  36  minute  first  dive,  5797  clicks  were  produced 
by  the  whale.  Approximately  97%  of  these  clicks  were 
foraging  clicks  with  an  Inter-Click  Interval  (ICI)  between 
0.15  and  1  sec  [3].  The  mean  foraging  click  ICI  was  0.31 
sec  (std=0.05).  This  is  in  agreement  with  Johnson  [4],  whom 
also  observed  a  regular  click  ICI  0.37  seconds  for  a  M. 
densirostris  in  the  Canary  Islands.  The  probability  density 
function  (PDF)  of  the  foraging  clicks  ICI  for  the  first  dive  of 
M.  densirostris  has  a  broad  peak  and  is  only  slightly 
asymmetric,  with  gradual  decay  at  higher  ICI’s  similar 
(Figure  8).  Zimmer  et  al.  [3]  noted  a  similar  asymmetric 
PDF  for  Ziphius  cavirostris ,  although  they  also  had  a  very 
sharp  peak  in  0.38  sec  ICI  rather  than  the  broad  peak 
between  0.29  and  0.32  seconds  observed  here. 


Figure  8:  Probability  density  function  for  foraging  click  ICIs 

The  FFT  detector,  implemented  with  a  noise  variable 
threshold  of  34,  was  able  to  detect  49%  of  the  clicks  on  at 
least  one  hydrophone.  The  matched  filter  detector, 
implemented  using  a  threshold  of  28.7  was  able  to  detect 
92%  of  the  clicks  on  at  least  one  hydrophone.  The  filter 
template  was  a  M.  densirostris  click  recorded  on  an  AUTEC 
hydrophone  from  a  previous  year.  On  each  hydrophone,  the 
matched  filter  detector  performed  significantly  better  than 
the  FFT  detector  (Figure  9). 


Figure  9:  Detection  effacy  per  hydrophone 

3.2  Localization 

While  the  majority  of  the  clicks  were  detected  on  at  least 
one  hydrophone  using  the  matched  filter,  three  dimensional 
localization  was  still  difficult  due  to  the  need  for  at  least  4 
TDOAs  between  the  DTag  and  hydrophone  as  input.  Of  the 
5767  clicks  associated  with  the  hydrophones,  only  1%  were 
detected  on  four  or  more  hydrophones.  More  commonly, 
the  clicks  were  detected  on  only  one  hydrophone  (44%), 
two  hydrophones  (36%),  or  three  hydrophones  (11%). 
Only  8%  of  the  clicks  were  not  detected  at  all.  As  a  result, 
the  TDOA  trends  were  interpolated  as  discussed  in  section 
2.6  prior  to  input  into  the  hyperbolic  multilateration 
algorithm.  In  addition,  depth  from  the  DTag  was  also  used 
as  an  initialization  parameter  to  provide  better  convergence 
of  the  solution. 

The  577  localizations  estimated  using  at  least  3  measured 
TDOAs  and  only  one  interpolated  value,  were  used  to 
“ground- truth”  the  track  estimated  by  the  DTag  alone  [1] 
(Figure  10).  Usually,  the  DTag  tag  on  and  tag  off  positions 
are  known  and  can  be  used  as  absolute  start  and  end 
positions.  However,  in  this  case  the  tagging  vessel  had  to 
leave  the  range  due  to  the  presence  of  range  operations  and 
the  tag  off  position  is  unknown.  In  order  to  “ground  truth” 
the  DTag  track,  small  user-chosen  sections  of  the  track  were 
individually  fitted  to  time- synchronized  3D  localizations  by 
adjusting  the  swim  speed  to  a  least-squares  match  (Figure 
ii). 
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Figure  10:  Original  DTag  Kalman-filtered  track  and  corrected 
DTag  track,  each  grid  square  is  2  km  x  2  km. 


Figure  12:  Swim  speed  estimated  from  localizations 


3.3  Detection  Range 


Figure  ll.DTag  track  corrected  based  on  3D  hyperbolic 
localizations,  tick  marks  at  200-m  increments 

Swim  speed  was  independently  estimated  using  the  x,y 
hyperbolic  localization  position  and  the  time  synchronized 
depth  from  the  DTag.  The  probability  density  function  of 
swim  speed  is  shown  in  Figure  12.  A  swim  speed  of  1.25 
m/s  was  observed  the  most,  approximately  12%.  The  swim 
speed  probability  peaked  rapidly  and  decayed  slowly  with 
increasing  swim  speed.  This  method  has  great  potential  for 
error  due  to  the  interpolation  of  between  1  and  3 
hydrophones  necessary  to  arrive  at  the  x,  y  location. 


The  3D  localizations  created  using  the  matched  filter 
detection  data  were  used  to  determine  the  range  from  the 
whale  to  the  hydrophone  for  each  detection  method. 
Clicks  from  each  hydrophone  determined  to  be  valid  on  the 
basis  of  their  TDOA  with  the  DTag  were  used  to  estimate 
range  and  bearing  hydrophone.  The  maximum  detection 
range  for  both  methods  was  approximately  6500  m, 
significantly  greater  than  previously  estimated  [10].  The 
whale  was  traveling  generally  in  a  north-east  direction,  but 
was  observed  to  turn  at  various  times  in  all  directions 
(Figure  13).  The  off-axis  aspect  angle  between  the  caudal- 
rostral  axis  of  the  tagged  whale  and  the  receiving 
hydrophone  was  determined  by  subtracting  the  bearing 
angle  from  the  whale  to  the  hydrophone  from  the  heading 
measured  by  the  DTag.  The  detection  range  as  a  function 
of  off-axis  aspect  angle  is  depicted  in  Figure  14  and  Figure 
15.  The  majority  of  the  clicks  detected  at  far  ranges  were 
within  ±30  degrees.  With  decreasing  range,  a  greater 
number  of  clicks  were  detected  further  off-axis.  While  a  -3 
dB  beam  width  of  6°  has  been  suggested  for  Z.  cavirostris 
by  Zimmer  [3],  M.  densirostris  may  be  less  directional  due 
to  their  smaller  body  size  and  potentially  smaller  source 
aperature  [10].  From  these  figures,  it  is  evident  that  the 
matched  filter  detector  significantly  outperforms  the  FFT 
detector  at  longer  ranges. 


90  f 


Figure  13:DTag  heading  for  dive  1,  Due  North  =  0,  the  radius  axis 
is  probability  density  (%) 


whale  to  hydrophone  indicates  the  ability  to  detect  signals 
significantly  off-axis  at  lesser  ranges  and  out  to  far  ranges 
when  close  to  the  axis.  For  Dive  1,  92%  of  the  clicks 
produced  by  the  tagged  whale  were  detected  at  least  one 
hydrophone  within  the  array.  This  combination  of  long 
detection  ranges  and  increased  probability  of  detection  with 
the  matched  filter  indicates  that  wide-baseline,  broadband 
arrays,  such  as  at  AUTEC,  provide  an  excellent  opportunity 
for  long-term  monitoring  of  beaked  whale  populations  and 
successful  passive  acoustic  based  mitigation. 
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